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An optical procedure in the context of continuous variables to verify bipartite entanglement with¬ 
out destroying both systems and their entanglement is proposed. To perform the nondestructive 
verification of entanglement, the method relies on beam-splitter and quantum nondemolition (QND) 
interactions of the signal modes with two ancillary probe modes. The probe modes are measured 
by homodyne detections, and the obtained information is used to feedforward modulation of signal 
modes, concluding the procedure. Characterizing the method by figures of merit used in QND pro¬ 
cesses, we can establish the conditions for an effectively quantum scheme. Based on such conditions, 
it is shown that the classical information acquired from the homodyne detections of probe modes 
is sufficient to verify the entanglement of the output signal modes. The processing impact due to 
added noise on the output entanglement is assessed in the case of Gaussian modes. 
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I. INTRODUCTION 

Entanglement is one of the most fundamental resources 
for performing processes in quantum information and 
computation. Besides the technological possibilities, the 
entanglement challenges our understanding of the quan¬ 
tum world and its connection with classical physics [1-3]. 
Recently, many experiments have accomplished quantum 
communication protocols sending light signals over dis¬ 
tances of hundreds of kilometers [4-6]. In these exper¬ 
iments, the entanglement was an essential part. It has 
also been studied as the use of entangled signals over 
long distances may increase the applicability of quantum 
cryptography protocols [7]. Thus it is very natural to 
devise stages along the transmission of entangled signals 
that verify if such signals are really entangled, without 
destroying or excessively disturbing them during the ver¬ 
ification processes. In other words, for future quantum 
communications, nondestructive certification protocols of 
entanglement will be required, ensuring the use of entan¬ 
glement for subsequent processes of quantum informa¬ 
tion. Studies of nondestructive entanglement verification 
or analysis have been done in the framework of discrete 
variable systems, such as single photons [8-10[. 

Differently from the previously cited studies, quantum 
communications using entangled signals may also be car¬ 
ried out in the context of continuous-variable systems, 
e. g., bright beams. The light beams may be regarded 
as oscillation modes, such that the states are vectors of 
an infinite-dimensional Hilbert space and observables are 
continuous spectrum operators, analogously to the posi¬ 
tion and momentum operators of the quantum harmonic 
oscillator [11]. For light beams, we consider amplitude 
and phase operators, also called quadratures. The re¬ 
search field of continuous-variable systems is very active 
and has extensive literature. Examples of quantum in¬ 
formation protocols performed with continuous-variable 
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light modes are quantum teleportation [12-15], cloning 
[16-18], and telecloning [19-21]. In such cases, the en¬ 
tanglement is an essential ingredient, therefore its conser¬ 
vation and its verification are primordial to further more 
complex applications. 

Thus a minimally invasive measurement method to ob¬ 
serve the entanglement is desirable. Although all quan¬ 
tum measurement entails a back-action effect, we can 
measure the signal, in order to preserve some of the 
properties of its original state. In particular, a quan¬ 
tum nondemolition (QND) measurement is able to mea¬ 
sure an observable without disturbing it, at the expense 
of a back-action disturbance on the conjugate observ¬ 
able [22-25]. In quantum optics, the QND measurements 
were initially performed by coupling the signal and probe 
modes in nonlinear optical media, such as Kerr media, 
optical fibers (third-order nonlinearity) [26-28], and in 
optical parametric amplifiers (second-order nonlinearity) 
[29, 30[. Other proposals relied on feedforward modula¬ 
tion of signal modes and off-line squeezed probe modes 
[31-34], in which it is not necessary to strongly pump a 
nonlinear medium in line to the signal modes. The com¬ 
bination of off-line squeezed probe modes, linear optics, 
homodyne detection, and feedforward loop has inspired 
many other optical operations, such as squeezing [33, 35- 
37], implementation of the one-way computation [38, 39], 
realization of third-order nonlinear operation [40], and 
other varieties of QND interactions [41, 42]. 

In this paper we propose a nondestructive method to 
verify continuous-variable bipartite entanglement, which 
uses QND and beam-splitter interactions between the sig¬ 
nal modes and two other probe modes. After the interac¬ 
tions, the probe modes are measured by homodyne detec¬ 
tions, so that the obtained photocurrents serve both to 
calculate the entanglement condition, so as to modulate 
the signal modes by electro-optic feedforward modula¬ 
tion, as has been implemented in noiseless optical am¬ 
plifiers [43-45]. The proposed scheme has the benefit of 
not mixing up the output signal modes to each other, 
which would change the global properties of entangle¬ 
ment [46, 47]. Another relevant result is that the ob- 
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tained entanglement condition is valid to the output sig¬ 
nals, ensuring the quantum correlation properties result¬ 
ing from the process. The cost of the procedure, manage¬ 
able by the scheme parameters, is the addition of excess 
uncorrelated noise in both signal modes. 

This paper is organized as follows. In Section II, we de¬ 
scribe the entanglement verification procedure, in which 
is shown the quadrature transformations of the signal and 
probe modes, in each stage of the scheme. In Section III, 
we characterize the procedure by the well-known QND 
measurement criteria. In Section IV, based on the re¬ 
sults of the previous section, we present how we compute 
a sufficient entanglement condition to the two output sig¬ 
nal modes, by detecting both probe modes. The effects 
of the noise addition on the signal entanglement are as¬ 
sessed in Section V, where the entanglement degradation 
is obtained in the case of Gaussian systems. Finally, we 
discuss the results and possibilities of this scheme in Sec¬ 
tion VI. 


II. ENTANGLEMENT VERIFICATION 



Figure 1. Schematic setup for the entanglement verification 
of modes 1 and 2. The probe modes are modes A and B. 
HD: homodyne detection for quadrature operators xa and 
Pb', BS: beam-splitter coupling with transmittances Ti and 
T 2 \ QND: quantum nondemolition coupling with gains Qi 
and Q 2 \ AM-PM: feedforward modulators. rrixA and rrip^ 
are the photocurrents obtained in detectors. These signals 
are used in both the feedforward process and to compute the 
entanglement condition. 


The goal of the process is to verify if a pair of sig¬ 
nal modes are entangled without destroying them. Since 
we consider the signal modes as continuous-variable sys¬ 
tems, we write their input quadrature operators as 
and pi for mode 1 and X 2 and p 2 for mode 2. In order to 
carry out the entanglement verification, two independent 
auxiliary beams must be introduced, characterized with 
input quadrature operators xa and pA for mode A, and 
xb and Pb for mode B. All operators obey the usual com¬ 
mutation relations, {xi,Pj\ = iSij^ for {i;j} = {1; 2; A; B} 
[11]. The auxiliary beams are used as probe modes, in¬ 
teracting with the signal modes and after measured by 
homodyne detections. Thus, the states of these probe 
modes A and B must be previously known. Following 
the theoretical proposals of previous articles [31-33], the 
probe modes must be prepared in strongly squeezed vac¬ 
uum states. In what follows, the squeezed quadratures of 
the probe modes are related to the operators pA and x b ■ 
First, each signal mode is coupled to each probe mode by 
ideal QND interactions. Interactions such as these have 
been performed by coupling beams in nonlinear optical 
media (see [24] and citations therein). However QND in¬ 
teractions were also performed using only linear optics, 
an auxiliary squeezed beam, homodyne detection, and 
feedforward modulation [33, 34] . As illustrated in Figure 
1, the beam pairs (1, A) and (2, B) are coupled by QND 
interactions, such that for modes 1 and A, with gain Gi: 


x[=xi, (1) 

'Pi=Pi-Gipa, ( 2 ) 

Xa = xa + GiXi, (3) 

p'a=Pa, (4) 


and for modes 2 and B, with gain G 2 : 


x'2= X2+G2 Xb^ (5) 

P2=P2, ( 6 ) 

x'b=xb, (7) 

p'b=Pb-G2P2- (8) 


After these first two QND interactions, the modes in¬ 
teract again, crossing probe modes with signal modes. 
Both interactions operate as beam splitters with trans¬ 
mittance Ti for modes 1 and B and transmittance T2 for 
modes 2 and A. Thus from equations (l)-(8), we have to 
modes 1 and B: 

x'l = \/t^xi - \/l- TiXb, ( 9 ) 

Pi = V^ipi - Gipa) - \/l - Ti{pb - G2P2), (10) 

x'b = V^XB + \/l - Tiii, (11) 

Pb = V^iiPB - G2P2) + v'l - Ti{pi - Gipa), (12) 
and to modes 2 and A: 

X2 = \/%{x2 + G2Xb)-Vi^^{xa + GiXi), (13) 
P2 = \/^P2 - a /1 - T2PA, (14) 

x'a = \/T^{xa + GiXi) + \/l — T2{X2 + G2Xb), (15) 
Pa = V^PA + \/l - T2P2- (16) 

These first two steps are necessary for the probe modes 
to obtain sufficient information from the signal modes. 
Then the probe modes are measured by homodyne detec¬ 
tion processes, providing classical signals (photocurrents) 
sufficient to compute and verify if the signal modes are 
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entangled. However, as we see in equations (9), (10), 
(13) and (14), the signal modes are affected by interac¬ 
tions. These perturbations can be corrected a posteriori 
with feedforward modulations, using phase and ampli¬ 
tude electro-optic modulators. Setting up the local oscil¬ 
lators of the homodyne detections, we select the quadra¬ 
tures x'^ and p'g to measure. With this choice, we obtain 
access to all signal quadrature operators, found in linear 
combinations in equations (12) and (15). In Section IV, 
we show such measures are sufficient for the entanglement 
verification. 

As the currently available detectors can achieve effi¬ 
ciencies above 99% [34], we will discard the noise from 
the detection process, so that we will focus only on the in¬ 
herent aspects of the procedure. However, the noise from 
detector imperfections can be calculated, which would 
add vacuum fluctuation terms in our derivations. 

The photocurrents generated in the detectors, rripg and 
nixA ) amplified with gains ki and k 2 for electro-optic 
modulators, so that the following signal quadratures are 
transformed as 


P'l = Pi + kirrip^ ( 17 ) 


III. QND CHARACTERIZATION 


As we are studying a procedure that must conserve 
some property of the signal modes, we must check it re¬ 
garding the features of a QND process. In early articles 
on QND measurement in optical systems, quantities were 
settled to characterize a device if it works as a noise¬ 
less amplifier and as a quantum state preparation (QSP) 
[24, 48-52[. To assess these features, we must consider 
quantities connecting statistical properties of the input 
and output modes, to both signal and probe pairs. The 
noise inserted in the system can be quantified by signal- 
to-noise ratios of the input signal, i?*", output signal, 
and output probe, The transfer coefficients 

from the input signal to the output signal are given by 


Ts 


R° 


m 


ym 

s 


-f Ns ’ 


(23) 


and from the input signal to the output probe by 


% = 


TDOUt 

Jlp 


yin 

* s 


Vr+A^p’ 


(24) 


and 

x'i'= X2 + k2mxA- (18) 

The respective conjugate quadratures remain unchanged. 
The gains ki and k 2 must be tuned in a way that the 
crossed terms between modes 1 and 2 are canceled. 

In conclusion, we can implement squeezing operations 
with gains Ti and T 2 onto signal modes 1 and 2 (not 
shown in Figure 1) , so that we obtain the output modes 


xT*=Xi-giXB, (19) 

pT^=Pi-Gipa, ( 20 ) 

x^^*" = X 2 + G 2 XB, ( 21 ) 

( 22 ) 


where gi = y^(l — Ti)/Ti and g 2 = y^(l — T 2 )/T 2 . The 
last squeezing operations onto modes 1 and 2 are not 
critical, because the entanglement is invariant under lo¬ 
cal linear unitary Bogoliubov operations [46, 47]. In 
equations (19)-(22), we maintain the terms with strongly 
squeezed quadratures, pA and xb, although their vari¬ 
ances tend to vanishing. At this point, it is interesting 
to present every possible resulting noise inherent to the 
scheme, reminding one that losses due to the efficien¬ 
cies of the detectors are not being considered. In fact, 
the variances of pA and xb are smaller the larger the 
squeezing in the probe modes. Nevertheless we need to 
know the scales of gi and Gi, i = {1,2}, before ruling 
out negligible terms. In the next section, the conditions 
to a genuine QND process will be studied, in which it 
is shown that the parameters gi and Gi can be tuned to 
optimize the scheme, so that some terms are eventually 
negligible. 


where V/” is the input signal quadrature variance, and 
Ns and Np are the equivalent input noises related to sig¬ 
nal and probe inputs, respectively. Another quantity of 
interest is the conditional variance of the output signal 
related to the measured output probe, given by 

WQSP = Vr-^-^, (25) 

where and 1}}“* are the signal and probe output 

quadrature variances, respectively, and Gs'^ is the sym¬ 
metrized covariance between former quadratures. Ac¬ 

cording to early articles [48-52], a fully QND process 
must simultaneously meet the following conditions: 

”7^ + 7); > 1, (26) 

indicating the noiseless amplifier property (quantum op¬ 
tical tapping), and 

Wqsp < 1, (27) 

indicating the quantum state preparation property. 

In the case of the entanglement verification, there are 
two output signals, each one with two quadratures, and 
two output probe quadratures. So we must calculate 
transfer coefficients (23) and (24) and conditional vari¬ 
ances (25) to a bipartite signal and a bipartite probe. 
That entails more combinations among the quadrature 
operators, implying more transfer coefficients and condi¬ 
tional variances to be considered. At first, we can seek for 
all combinations of quadrature operators between signal 
and probe systems. On the other hand, transfer coef¬ 
ficients crossing quadratures do not exist (for example, 
there is no Np^^ related to Xi). Moreover, a direct veri¬ 
fication unveils that conditions (26) and (27) cannot be 
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satisfied in all existing combinations of quadrature oper¬ 
ators. However, to meet a quantum regime, it is sufficient 
to comply only with conditions (26) and (27) related with 
the signal quadratures preserved in the QND procedure. 
Therefore, all these considerations restrict the relevant 
transfer coefficients and conditional variances. For ex¬ 
ample, we can choose the signal quadrature operators 
xi and p 2 to be preserved. Thus the equivalent noises, 
introduced in systems, are 

(28) 

Np,=gl{{^PAf), (29) 

= ^((Apb)^) + ((Api)^) 

- |L(i{Api,Ap2})+ (^)'((Ap^)2), (30) 

= ^((Axa)^) + ((Ai^)^) 

+ ^(i{Aii, Ai2}) + ((Axs)^), (31) 

where = Oi — {(Di) and {Oi) = Tv{Oip)^ such that Ot 
is some operator distinguished by index i, and p is the 
density matrix of the whole system. With expressions 
(28)-(31), we can find the respective conditions (26). Af¬ 
ter simple algebra, such conditions are rewritten as 


and 


< {{Ax,)y 


< {iAp2)Y 


(32) 


(33) 


So it is clear that sufficiently small values of gi and g 2 
and sufficiently large values of Gi and G 2 will achieve 
a quantum optical tapping regime. The ideal situation 
is obtained when 51 ; 52 —^ 0 and G'i;G '2 —>■ c», in which 
perfectly noiseless amplifiers are achieved. 

Conditional variance (25) applied to output modes un¬ 
folds in two other quantities: 


Wqsp{xi,xa) = (( Axr ‘)^) - 


|(i{AirSAx^“‘})|= 


((Ax7*)2) 


and 


(34) 


Wqsp{p2,Pb) = {{Ap777 - 


|(^{Apr.Apg-*})p 

((Apr)") 


(35) 

It is possible to find conditions to the quantum state 
preparation regime with both expressions (34) and (35) 
simultaneously. A limit case can be obtained, consid¬ 
ering strongly squeezed input probe modes, such that 
{{Axb7)', {{^Pa)7 0, and based on previous consid¬ 

erations, taking 51152 -A 0, we notice that the quantum 
state preparation is attainable if, from expression (34), 


Gl>il- 


1 


((Aii) 2 ) 


((Ax^)"). 


(36) 


As (1 — l/((Aii)^)) < 1 , to any physical value of 
((Axi)^), a stricter inequality is more interesting, 

Gl > {{Axa)7- (37) 

Similarly, from expression (35), we can obtain another 
inequality, 

Gl > ((Aps)"). (38) 

Both inequalities (37) and (38) can be fulfilled simulta¬ 
neously, therefore quantum state preparation regimes are 
feasible for reasonable values of the parameters of the op¬ 
tical device, independently of the input beam properties. 

Considering other combinations of signal quadrature 
operators to be preserved in the QND procedure, we can 
seek other conditions to 51 , 52 , Gi, and G 2 , analogously 
to previous analysis. Such cases are very similar and do 
not add new information. For the example studied, we 
can safely approximate the output modes to 


= xi, (39) 

pT^=Pi-Gipa, (40) 

= ^2 + G 2 XB, (41) 

Pr=P 2 . (42) 


We can notice in the output signals that each mode has a 
preserved quadrature, featuring a QND process. On the 
other hand, each mode has a conjugate quadrature added 
by terms from probe modes. As the input probe modes 
are independent, these terms produce phase-sensitive un¬ 
correlated noise, inevitably disturbing the signal modes. 


IV. CALCULATING THE ENTANGLEMENT 
CONDITION 

Besides feedforward modulation, the photocurrents are 
also used to calculate the entanglement condition of the 
signal modes. Duan et al. [47] have found a sufficient 
entanglement condition in continuous-variable systems, 
based on EPR-like operators. Later, other works have ex¬ 
tended this condition for more general operators [53, 54[ . 
Following these authors, consider operator combinations 
such as 

11 = 01 X 1 - 1 - 02 X 2 , (43) 

V = bipi + 62 P 2 , (44) 

where [xi,pj] = i5ij, (i]j = 1 ; 2 ), and Oi, 02 , 61 , and 
&2 are arbitrary constants. It is possible to show that 
a sufficient condition of continuous-variable bipartite en¬ 
tanglement is 

((All)2) + ((AD)2)<|ai6i| + 10262 ], (45) 

valid to Gaussian or non-Gaussian systems. 

On the other hand, the variances of the photocurrents 
generated in the homodyne detection can be written as 

{{Amx^f) = KaT2 [((Ao°“*)2) + ((Ax^)^)] , (46) 
((Ampj2) = KbT, [((Ah°“*)2) + ((Apb)2)] , (47) 
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where + 92xT* and {)°“* = - G 2 PT* 

are defined from equations (39) to (42). Ka and Kb are 
factors that depend on the overall detector efficiencies 
and the conversion circuitry of the photocurrents, such 
that these factors can be related by Ka = ( 52 /^ 2 )^ and 
Kb = {gi/kif‘. Except for ((A'u°“*)^) and ((Af)°“‘)^), 
all other terms of equations (46) and (47) are measured 
or previously known. So both ((Au°“‘)^) and ((A{}°“*)^) 
can be calculated and compared with equation (45), iden¬ 
tifying the arbitrary constants with the parameters of 
the apparatus, i. e., ai = Gi, 02 = 321 ^1 = 3i) and 
&2 = — G 2 . Therefore the entanglement condition (45) 
can be calculated from the detected probe modes and 
from parameters of the scheme. We can also notice that, 
to sufficiently squeezed probe modes, the photocurrent 
variances are directly proportional to the EPR-like oper¬ 
ator variances, i. e., ((Amxyi)^) — KaT 2 {{Au°'^*Y) and 
{{Ampg)'^) ~ ifTsTi((Af)““*)^), so that the photocurrent 
measurements provide a direct way to verify the entan¬ 
glement. 

An important result of this method is that the condi¬ 
tion of entanglement (45), calculated with expressions 
(46) and (47), is exactly valid for the output signal 
modes, namely, we can certify that the signals resulting 
from the scheme are entangled, regardless of limitations 
or scheme losses. From a practical point of view, we are 
receiving two signal modes to verify the quantum proper¬ 
ties of their correlations, so that we can nondestructively 
maintain them for use in other processes, and still be able 
to repeat it. Inevitably the scheme has a cost, which is 
the addition of phase-sensible noise, degrading the signal 
modes. The effects of this degradation on the entangle¬ 
ment are discussed in the next section. In an idealized 
situation, the probe modes would have a squeezing pa¬ 
rameter tending to infinity, so the signal modes could 
have their entanglement checked without any degrada¬ 
tion, perfectly preserving each quadrature of the signals 
as well. 


V. ENTANGLEMENT DEGRADATION 


To assess the effects that the presented scheme cause on 
the signal mode entanglement, we will restrict this anal¬ 
ysis to the case of Gaussian beams [55]. The continuous- 
variable systems restricted to Gaussian states are fully 
described by the first statistical moments of the dynam¬ 
ical operators, i. e., Oi = {Oi), and by the second sta¬ 
tistical moments, that can be arranged in a covariance 
matrix, M, whose entries are M^j = i({A(!)j, AOj}). 
With a suitable choice of quadrature basis, we can write 
the input covariance matrix as 


Min 


/ ni 0 c 0 \ 

0 ni 0 fc I 

C 0 712 0 I 

Vo k 0 712 / 


(48) 


After all procedures, the covariance matrix of the signal 
modes becomes 


Mout 


/ ni 0 c 0 \ 

0 Til -|- di 0 k \ 

C 0 712 + ^2 0 I ’ 

Vo k 0 772 / 


(49) 


in which we can observe the presence of uncorrelated ex¬ 
cess noises di = G\{{ApaY)) and ^2 = G\{{AxbY)), 
generated from probe operators in equations (40) and 
(41). These noises spoil the input entanglement. This 
can be seen by calculating the logarithmic negativity [56- 
58]: 

£'Ar(p) = max[0, —Lnf/-], (50) 

where is the smallest symplectic eigenvalue of the 
partially transposed bipartite state. This quantity can 
be calculated from symplectic invariants of the partially 
transposed system: 

/A — \/ A^ — 4 det M 

= Y-: (51) 

where A = 771(771 -|- di) + 772(772 -I- £^ 2 ) — 2cfc and det M = 
[ 771(772 -I- d 2 ) — c^][772(77i -|- di) — fc^], if we use matrix 
(49). Hence we find the expression very complex. To 
simplify the problem, we consider input signal modes as 
two-mode squeezed vacuum, so 77i = 772 = cosh(2?’) and 
c = —k = sinh(2?’). With these substitutions, we plot 
Figure 2. 



Figure 2. Plot of logarithmic negativity to a two-mode Gaus¬ 
sian state characterized by matrix (49) as the function of 
excess noises di and d 2 . The input modes are two-mode 
squeezed pure states, with squeezing parameter r — 1. 


One can see in Figure 2 that, for any nonzero values 
of di and d 2 , the logarithmic negativity is monotonically 
decreasing. Such effect is the verification process cost. 
The optimization is reached minimizing the variances of 
Pa and Xb, i- e., increasing the squeezing of the probe 
modes. In addition, the gains Gi and G 2 must be limited, 
obeying conditions (32), (33), (37), and (38). 
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VI. DISCUSSION 

In this paper, a nondestructive scheme for bipartite en¬ 
tanglement verification in the framework of continuous- 
variable systems was shown. To such task, a suitable 
choice of QND and beam-splitter interactions between 
the pair of signal modes and a pair of probe modes is 
necessary, followed by measurements of the probe modes 
and feedforward modulations of the signal modes. All 
of these processes are feasible with current technologies, 
therefore it can be performed in experimental demonstra¬ 
tions or implemented as a built-in step in a larger com¬ 
munication protocol, in which it is necessary to certify 
that two signals are entangled, while they are used in an¬ 
other further step. Some studies have already been done 
with similar purposes in the context of discrete variable 
systems [8-10]. So this paper fills a gap for continuous- 
variable systems. 

This method is based on a sufficient entanglement con¬ 
dition, therefore some entangled states cannot be de¬ 
tected. However, entangled states that do not satisfy 
inequality (45) are fragile when subjected to Gaussian 
attenuation, as already shown in the articles by Barbosa 
et al. [59, 60]. So these fragile states would not be in¬ 
teresting for long-distance implementations. Moreover, 
maximally entangled states or near them are required 
for many quantum information protocols. Such states 
are addressed by the scheme. 

According to the criteria of QND measurement char¬ 
acterization [48—52] , we assess the quantum properties of 
the entanglement verification scheme. One notice that 
there are many possibilities to calculate the transfer co¬ 
efficients and the conditional variances, using different 
quadratures of the bipartite modes. In this paper, we 
calculate all relevant quantities to certify the QND prop¬ 
erties, although we do not present a systematic method 
for finding them. Therefore a QND multipartite device 
characterization would be a relevant theoretical develop¬ 
ment for further research. 


The entanglement of the signal modes can be checked 
by the presented method, but it cannot quantify the en¬ 
tanglement. Such limitation exists because the photocur¬ 
rents, measured by homodyne detectors, provide infor¬ 
mation only about the variances of the EPR-like opera¬ 
tors, aiXi -|- 02^2 and bipi + & 2 P 2 - That is insufficient to 
have a measure of entanglement, e. g., logarithmic neg¬ 
ativity, although it is sufficient to detect entanglement. 
However, new strategies may lead to quantifying the en¬ 
tanglement. We can expect that more complex signal- 
probe interaction configurations reach these goals. Un¬ 
like squeezed vacuum modes, using non-Gaussian modes 
as probe modes could also lead to more promising results, 
as already noted in articles concerned with the trade-off 
between information and disturbance caused by measure¬ 
ments in continuous-variable systems [61, 62]. All these 
considerations show many future research possibilities. 

We can notice that the entanglement verification is 
deeply connected with the eavesdropping in quantum 
cryptography [63]. While an eavesdropper wants to make 
a necessarily imperfect copy of the signal sent between 
two communication stations, in the proposed entangle¬ 
ment verification, it looks to observe the correlation be¬ 
tween the two signal modes, without necessarily copying 
and measuring the physical states. The similarity be¬ 
tween the two processes is that both are extracting infor¬ 
mation from the signal, and both are adding unavoidable 
perturbation, in the case of this paper, adding uncorre¬ 
lated phase-sensible noise. In both processes, the mini¬ 
mization of the perturbation is the condition to its op¬ 
timal accomplishment. Further studies may be devoted 
to the details and explanations of the relationship be¬ 
tween the QND verification of multipartite correlations 
and cloning and eavesdropping. 
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